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ABSTRACT 

High-resolution time series of sunspots have been obtained with the Swedish lm Solar Telescope 
between 2003 and 2006 at different locations on the solar disc. Proper motions in seven different 
active regions have been studied. The analysis has been done by applying local correlation tracking 
to every series of sunspots, each of them more than 40 minutes long. The sunspots' shapes include a 
different variety of penumbral configurations. We report on a systematic behaviour of the large-scale 
outflows surrounding the sunspots, commonly known as moat flows, that are essentially present only 
when preceded by a penumbra not tangential but perpendicular to the sunspot border. We present 
one case for which this rule appears not to be confirmed. We speculate that the magnetic neutral 
line, which is located in the vicinity of the anomalous region, might be responsible for blocking the 
outflow. These new results confirm the systematic and strong relation between the moat flows and 
the existence of penumbrae. A comparative statistical study between moats and standard granulation 
is also performed. 

Subject headings: Sun: sunspots - Sun: granulation - Sun: photosphere 



1. INTRODUCTION 

The Sun's strong magnetic behaviour is clearly re- 
vealed through its most prominent visible manifestation, 
the sunspots. The newest generation of solar telescopes 
and the lastest restoration techniques have greatly in- 
creased the spatial resolution of the images and therefore 
the variety of details and tiny structures in and around 
sunspots that can be resolved. Nevertheless, it is still not 
well understood how these structures are formed, evolve 
and affect the photosphere surrounding them. Many ap- 
proaches have been done in this direction such as the 
identification of small magnetic elements called "Moving 
Magnetic Features" (MMF) travelling radially outwards 
while immersed in an anullar cell around the sunspot of 
strong radial outflows known as a "moat" (Sheeley 1972; 
Harvey & Harvey 1973). For a recent overview of the 
literature on MMFs see the introduction of Hagenaar & 
Shine (2005). The sunspot moat thus defines an orga- 
nized horizontal flow pattern which ends quite abruptly 
at a distance that can be comparable to supergranular 
sizes or even larger. Sobotka & Roudicr (2007) study the 
characteristics and temporal evolution of moats at two 
heights in the atmosphere for a large sample of sunspots 
with different sizes, shapes and evolutive stage. A simi- 
lar work studying the properties of moats was previously 
performed by Brickhouse & Labonte (1988). The aver- 
aged moat flow velocity reported ranges from ^0.4 to ~1 
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kms . 

New findings (Sainz Dalda & Martinez Pillet 2005) 
have shown that the penumbral filaments extend be- 
yond the sunspot boundary entering the region domi- 
nated by the moat flow where the MMF activity is de- 
tected. Thus, the temporal average of magnetograms 
has unveiled the existence of moat filaments: horizon- 
tal, filamentary structures coming from the penumbra 
and reaching the photosperic network as an extension of 
penumbral filaments. Moreover, some MMFs have been 
found (Sainz Dalda & Martinez Pillet 2005; Ravindra 
2006) starting just inside the sunspot boundary in its 
way out from the sunspot. Kubo ct al. (2007) have es- 
tablished a relationship between the vertical components 
(spines) of the magnetic field in the so-called uncombed 
structure of the penumbra and MMFs observed in moat 
regions. 

In a recent letter, Vargas Domfnguez et al. (2007) stud- 
ied the moat flow in a complex active region. They found 
that no outflow was detected in the granulation next to 
umbral cores that lack penumbrae. Outflows were only 
found in the granulation regions adjacent to the penum- 
brae in the direction following the penumbral filaments. 
Granulation regions located next to penumbral sides par- 
allel to the direction of the filaments show no moat flow. 

The aim of the present work is to extend the study 
by Vargas Domfnguez et al. (2007) to a larger sample of 
cases, that is to establish whether this moat-penumbra 
relation is systematically found in other active regions 
and how the granular convective pattern surrounding 
sunspots behaves. In doing so, we have used several 
high quality time series observed at high spatial resolu- 
tion. A total of 7 different sunspot series have been pro- 
cessed and analyzed. The sample includes sunspots with 
different penumbral configurations, varying from well- 
developed penumbrae to rudimentary penumbral mor- 
phologies. The paper describes first in § 2 how the ob- 
servations were performed and the restoration technique 
used to correct for atmospheric aberrations. Once the 
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time series are ready to analyze, § 3 presents some cal- 
culations and statistics of proper motions of structures 
surrounding the sample sunspots. We finally discuss the 
results and future work in § 4. 

2. OBSERVATIONS AND DATA PROCESSING 

The observations were obtained with the 
Swedish 1-m Solar Telescope (SST, Scharmer et al. 
2003a) on La Palma between 2003 and 2006. All 
observations benefited from the use of the SST adaptive 
optics system (Scharmer et al. 2003b) that minimized 
the degrading effects of seeing. Image post-processing 
techniques were applied to increase the homogeneity 
in the quality of the time series and to enhance the 
image quality over the whole field-of-view (FOV). For 
the 2003 data set we applied Multi-Frame Blind Decon- 
volution (MFBD) using the implementation developed 
by Lofdahl (2002). The MFBD code was succeeded 
by the Multi-Object Multi-Frame Blind Deconvolution 
code (MOMFBD, Van Noort, Rouppe van der Voort 
& Lofdahl (2005)) which employs multiple objects and 
phase-diversity. MOMFBD was applied to the data 
sets after 2003. For all the time series we present in 
this paper, the seeing conditions were generally very 
good and sometimes excellent. A large fraction of the 
restored images in the different time series approach the 
diffraction limit of the telescope. 

After image restoration, we applied standard tech- 
niques to the time scries including correction for the di- 
urnal field rotation, rigid alignment and de-stretching to 
correct for seeing-induced image warping. A p-modes 
filter was also applied to the series (threshold phase ve- 
locity 4 kms' 1 ). 

Table 1 presents details for the different active region 
targets, Table 2 gives some details on the restored time 
series. Below we provide more detailed information on 
the different data sets. 

2.1. SI: AR440, 22-Aug-2003 

An interference filter centered at the continuum at 

436.4 nm (FWHM 1.1 nm) was used. In the remainder, 
this filter will be referred to as "G-cont" filter. Every 
25 s, the 3 highest contrast images were selected and 
stored to disk. These 3 images were used for MFBD 
processing. This sunspot is shown in Figure 2. 

2.2. S2: AR608, 10-May-2004 

An interference filter centered at the G-band at 

430.5 nm (FWHM 1.2 nm) was used. In the remainder, 
this filter will be referred to as "G-band" filter. This ob- 
servation involved 2 cameras set-up as a phase-diversity 
pair: one in focus and one slightly out of focus. 20 expo- 
sures from each camera, acquired over the course of 7.5 
seconds, were used for MOMFBD processing. During 
processing of the time scries, a few bad quality images, 
suffering from too much blurring, were dropped. The 
exposure time was 11 ms. This sunspot is shown in Fig- 
ure 3. 

2.3. S3: AR662, 20-Aug-2004 

This time series was observed with the G-cont filter. 
The MOMFBD processing was done using 10 exposures 
(acquired during 7s) from a phase-diversity pair of cam- 
eras. For the central part of the FOV, centered on the 



sunspot, the MOMFBD processing was extended to in- 
volve a G-band phase-diversity pair of cameras. After 
MOMFBD restoration, the best quality image from 3 
subsequent images was selected for further processing of 
the time series. The exposure time was 11.3 ms. This 
sunspot is shown in Figure 9. 

2.4. S4: AR662, 21-Aug-2004 

This time series was processed in a similar manner as 
S3. After MOMFBD restorations, the worst 14% of the 
images were dropped from the time series. This sunspot 
is shown in Figure 4. 

2.5. S5: AR789, 13-Jul-2005 

For this time series, MFBD processing was applied to 
18 subsequent G-band images, acquired during 10 s. Im- 
ages that were too blurred and when the AO system 
was not actively compensating, were dropped from the 
MFBD sets. This sunspot is shown in Figure 5. 

2.6. S6: AR813, 4~Oct-2005 

This time series was observed with an interference fil- 
ter centered on the spectral region between Ca K and 
Ca H at 395.4 nm (FWHM 1.0 nm). For the MOMFBD 
processing, 20 exposures (acquired during 8.3 s) from a 
phase-diversity pair of cameras were used. In addition, 
simultaneous exposures from two other cameras were 
used for the restoration: one equipped with a narrow- 
band filter centered in the Ca H wing at 396.5 nm and 
one equipped with a narrow-band filter centered on the 
Ca H core (396.8 nm). The exposure time was 11 ms. 
This sunspot is shown in Figure 6. 

2.7. S7: AR893, 10-Jun-2006 

This time series was observed with an interference filter 
centered at 630.2 nm (FWHM 0.8 nm). This is the pre- 
filtcr for the Lockheed SOUP filter. MOMFBD process- 
ing was applied to 400 exposures obtained during 19 s. 
The restorations included images from a phase-diversity 
pair of cameras with the wide band filter, and narrow 
band SOUP exposures in the Fe I 630.2 nm spectral line. 
This sunspot is shown in Figure 7. 

3. DATA ANALYSIS 

The high quality, stability and long duration of the re- 
stored sunspot series, enable us to follow the dynamics 
of the plasma around the sunspots. We have computed 
proper-motion velocity fields (horizontal velocities) em- 
ploying the local correlation tracking (LCT) technique 
(November & Simon (1988) using the implementation of 
Molowny-Horas & Yi (1994)). A Gaussian tracking win- 
dow of FWHM l'.'O suitable for tracking mainly granules 
has been used for all the series. The method produces a 
sequence of map-pairs each one describing, within a short 
time interval (tens of milliseconds) , the horizontal veloc- 
ity components — x and — y along the FOV. These maps 
of components have been averaged over 5 and 10 min pe- 
riods, and also over the respective total duration of every 
sunspot series (more than 40 min in all cases). From the 
averaged maps of velocity components, the distribution 
of velocity magnitudes over the FOV (flow map) has been 
easily derived. 

It is well documented in the literature that LCT- 
techniques in general produce some systematic errors 
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TABLE 1 

Sunspots Sample 



Name 


Active region 
NOAA 


Observing date 


Heliographic coordinates 




SI 


10440 


2003, Aug 18 


(S10, Wl) 


0.96 


S2 


10608 


2004, May 10 


(S5, E57) 


0.55 


S3 


10662 


2004, Aug 20 


(N13, E7) 


0.99 


Si 


10662 


2004, Aug 21 


(Nil, W5) 


0.99 


S5 


10789 


2005, Jul 13 


(N17, W32) 


0.88 


S6 


10813 


2005, Oct 04 


(S7, E37) 


0.78 


S7 


10893 


2006, Jun 10 


(Nl, E17) 


0.95 



TABLE 2 
Restored sunspots time series 



Name 


N. images 


Cadence (s) 


Duration (min:s) 


SI 


128 


24.7 


52:47 


S2 


376 


7.5 


47:10 


S3 


144 


19.8 


47:20 


S4 


344 


8.0 


45:57 


S5 


240 


10.1 


40:12 


S6 


556 


8.7 


80:3 


S7 


124 


19.7 


40:26 



in the determination of displacements (see Figure 13 in 
November & Simon (1988)). Their significance depends 
on several factors like the width of the tracking window 
in relation to the size of the elements used as motion trac- 
ers, and the presence of large intensity gradients in the 
images. Also the nature of the interpolation algorithm 
to fix the position of the correlation maximum with sub- 
pixel accuracy may play a role. The LCT-method typi- 
cally underestimate the velocity fields. Tests performed 
by the authors of the LCT-code we are employing in this 
work lead to conclude that in the worst cases this under- 
estimation can amount to 20-30 % (Yi 1992; Molowny- 
Horas 1994). Even so, this drawback does not change 
the main conclusions of the present paper since we are 
not interested in fixing absolute values of velocities but 
rather in the detection of large-scale regularly organized 
flows around the sunspots in comparison with those typ- 
ically found in a normal granulation field. In both cases 
we use the solar granules as the tracers of motion and 
the same size for the tracking window so that in case of 
some bias it will affect similarly to both velocity fields. 

The sample of sunspots studied in this paper includes a 
variety of heliocentric positions on the solar disc as shown 
in Table 1 (with /i = cos 8, where 9 is the heliocentric an- 
gle). Away from solar disc center, the measured proper 
motions are in fact projections of the real horizontal ve- 
locities in the sunspot plane onto the plane perpendicu- 
lar to the line-of-sight (LOS). Thus, to evaluate the real 
horizontal velocities we have de-projected the velocities 
obtained from the LCT technique to a plane tangent to 
the solar surface at the center of the sunspot. 

In Figure 1 we consider the orthogonal coordinate sys- 
tem SX,SY,SZ (the Sunspot System, SS) with the SZ- 
axis perpendicular to the plane tangent to the solar sur- 
face at the sunspot location and the SX-axis tangent to 
the meridian of the solar sphere at this location. The 
figure also shows the coordinate system X,Y,Z (the Ob- 
serving System, OS) with a common origin, the Z-axis 
coinciding with the LOS and the X-axis pointing to the 
axis crossing the sun center in a direction parallel to the 
LOS. Note that the axis, X, SX, Z, and SZ are coplanar 
and SY and Y are collinear. The real horizontal veloci- 




FlG. 1. — Sketch showing two orthogonal coordinate systems for 
the projection analysis of the horizontal velocities: 1) The Sunspot 
System (SX,SY,SZ) with the Z-axis perpendicular to the solar sur- 
face at the sunspot location, and 2) the Observing System (X,Y,Z), 
where the observations are described, with the Z-axis coinciding 
with the line-of-sight. 

ties, v(v, (j)), are contained in the (SX, SY) plane whereas 
we observe their projections v'(v' , </>'), in the (X,Y) plane. 
Calculations based on Figure 1 lead to the following re- 
lationships (equations 1 and 2) between the magnitudes 
v and v' , and the azimuths <f> and 0', that will allow us 
to deproject our observations. 



= v 2 (sin 2 4> + 



cos 



cos 



tan ( 



tan i 



cos ( 



(1) 



(2) 



Since the X-axis and Y-axis do not necessarily coin- 
cide with our natural reference system (i.e the edges of 
the CCD frame), the measurement of the projected az- 
imuths, <j)' , requires the knowledge of the orientation of 
our FOV with respect to the X-axis that points to the 
solar disc center. 
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3.1. Flowmaps for different penumbral configurations 

As mentioned above, for every sunspot series, the map 
of horizontal velocities averaged over the whole time se- 
ries is evaluated. In order to coherently detect the moats 
and compare the statistics (§ 3.3) in all active regions of 
our sample, we have de-projected the observed velocity 
vectors onto the sunspot plane as described above. Over- 
laying the observed images, we construct maps show- 
ing in the granulation field only those projected velocity 
vectors having de-projected magnitudes above a certain 
threshold. In these maps, extensive organized outflows 
coming out from the sunspots can easily be identified 
so that we can outline masks delimiting in a qualitative 
way the area of the moats for the purpose of statisti- 
cal calculations and graphical representation. We have 
fixed the mentioned threshold to <~ 0.3 kms -1 as the 
best compromise to define the limits of the moats. This 
value should not be understood as an absolute (univer- 
sal) velocity threshold to define moats in general. As we 
have mentioned before the particular method used here 
to measure the displacements may underestimate the ve- 
locity magnitudes. Moreover, outside the moat one can 
easily find velocities larger than that. The point is that 
this particular threshold makes in our case more evident 
the existence of a well organized radial outflows around 
the sunspots in comparison with the rest of the FOV. 
Lower values (even zero) for the threshold do not extend 
the areas of organized flows around the sunspots but pro- 
duce maps with a very dense and noisy (exploding gran- 
ules everywhere) representation of arrows where the out- 
line of the frontiers of the moats becomes more difficult. 
Figures 2 to 6 show within the moats, only those velocity 
vectors with de-projected magnitudes above 0.3 kms -1 . 
In Figure 7 we extend this representation to the entire 
granulation field showing that large velocities are also 
present outside the moat. These velocities are generally 
grouped and associated with exploding granules. Our 
sunspot sample includes different penumbral configura- 
tions as a key factor to establish the moat-penumbra re- 
lation in a robust way. 

Firstly we focus on granulation regions that display 
moat flows. Close inspection of Figures 2 - 7 , and of 
the upper panel of Figure 9 reveals that the velocity vec- 
tors in the moats are oriented following the direction of 
the penumbral filaments. We observe that in all cases 
the moat flow direction lines up with penumbral filaments 
that are oriented radially with respect to the sunspot cen- 
ter. Nevertheless, the completeness of the velocity vec- 
tors (density of arrows) in the flow maps depends on 
the threshold previously imposed, as expected. This can 
be seen in Figure 5 where the left-hand side penumbra 
docs not seem to be strictly associated with large flows 
in the granulation region. However, this penumbral part 
has in fact an associated moat flow similar to the other 
penumbral regions which is not visible in the representa- 
tion since the magnitude of the velocities is slightly lower 
than the threshold of 0.3 kins" 1 . 

Next we consider granulation regions close to the 
sunspots that do not display moat flows. All sunspots 
of our sample, except S3, have parts of umbral core in 
direct contact with granulation regions without interven- 
ing penumbrae. In all these granulation regions we do 
not detect systematic large-scale outflows that fulfill our 



criteria for moat flows. This supports the conclusion that 
umbral core boundaries with no penumbra do not display 
moat flows. 

Finally, there are also granulation regions in the vicin- 
ity of penumbrae that lack significant moat flows. To 
study these cases in more detail we mark in Figures 2, 
4, 5 and 6 peculiar regions of interest with rectangular 
white boxes. These regions are represented as close-ups 
in Figure 8. In the three upper panels, the penumbral 
filaments display significant curvature such that they are 
not radially oriented with respect to the sunspot centre 
but they extend in a direction tangential to the sunspot 
border (as marked with black lines in the overview im- 
ages Figures 2, 4 and 6). The lower panel in Figure 8 
shows a penumbra extending radially from the umbra. 
The surrounding photosphere exhibits moat flows only 
in the direction of the penumbral filaments (also sec Fig- 
ure 5). These four examples lead to the conclusion that 
moat flows are not found in directions transverse to the 
penumbral filaments. 

In a few cases we found small pores in the vicinity 
of penumbrae in a region where we would expect moat 
flows. Such cases are seen around coordinates (19,29) 
in Figure 2 and coordinates (29,11) in Figure 6. For 
these regions the measured proper motions are less reli- 
able since we have a FWHM=l'/0 tracking window act- 
ing on a region of only a few arc seconds (~ 2-3) between 
the pore and penumbra. Another possibility is that the 
pores themselves are somehow blocking the large out- 
flows changing the expected behavior. In the case pre- 
sented in Figure 6 the most plausible explanation for the 
absence of moats is the presence of a neutral line close 
to coordinates (29,11) as we shall describe in § 3.2. 

The findings described above suggest a link between 
the moat flows and the Evershed flows in penumbrae. 
We come back to this relation in § 4. 

3.2. Neutral lines affecting the flows behavior 

Only one of the sunspots we study displays a complete 
regular and well-developed penumbra surrounding com- 
pletely the umbral core. Figure 9 shows the flow map 
calculated for this active region, plotting the horizon- 
tal velocities surrounding the sunspot with deprojected 
magnitudes > 0.3 kms" 1 . Large outflows arc not found 
in a certain part of the right-hand side of the spot in the 
black squared region in the upper panel of Figure 9. 

Following the findings of § 3.1, we would expect to find 
a moat flow in this region. The penumbral filaments arc 
oriented radially from the umbral core which suggest the 
presence of a moat flow in the granulation region in the 
direct vicinity. Even with a lower velocity threshold, no 
moat flow can be discerned in this region. Neverthe- 
less, when comparing with the magnetogram (see lower 
panel of Figure 9), we found an inversion in magnetic 
polarity just outside the right border of the sunspot: the 
magnetogram displays positive polarity (in white) for the 
sunspot but negative polarity for the small magnetic el- 
ements and pore just outside the penumbra. The rever- 
sal of polarity (or neutral line) is confined to a narrow 
region that roughly coincides with the sunspot border. 
The absence of large outflows following this penumbra is 
suggestively related to the presence of this neutral line 
which might somehow be acting as a blocking agent for 
the moat flow in this region. The position of the neutral 
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Fig. 2. — Sunspot SI: Map of the horizontal velocities inside the 
moat with de-projected magnitudes > 0.3 kms -1 . The sunspot 
exhibits a well developed penumbra in the lower region. The moat 
flow is clearly found in the same region following the penumbra fila- 
mentary direction and in no other region surrounding the sunspot. 
The coordinates are expressed in arc sec. The black bar at (0,0) 
represents 1.5 kms -1 for the projected velocities in all the flow 
maps hereafter. A peculiar region is white squared. 




Fig. 3. — Sunspot S2: Map of the horizontal velocities inside the 
moat with de-projected magnitudes > 0.3 kms -1 . This sunspot 
shows different penumbrae shapes all around it, from short to long 
and from wide to narrow penumbrae. Except for the places where 
the umbra is adjacent to the surrounding photosphere, we identi- 
fied moat flows all around the sunspot in roughly the same con- 
figuration of the penumbra. The coordinates are expressed in arc 
sec. 

line measured by LOS magnctograms is generally influ- 
enced by the location of the sunspot (8 angle). In this 
case since the sunspot is very close to the disc center we 
can claim a reliable determination of the neutral line. 

For the sunspot in Figure 6 there is also a neutral line 
crossing all along the right-hand border of the sunspot 




Fig. 4. — Sunspot S4: Map of the horizontal velocities inside 
the moat with de-projected magnitudes > 0.3 kms -1 . The figure 
corresponds to a sunspot with two distinct penumbrae. One tan- 
gential to the umbral core and the other one coming out radially 
from the right part of the umbra. The moat flow continues as pro- 
longations of the last one. The coordinates are expressed in arc 
sec. A peculiar region is white squared. 
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Fig. 5. — Sunspot S5: Map of the horizontal velocities inside the 
moat with de-projected magnitudes > 0.3 kms -1 . This sunspot 
presents very distinct areas with and without penumbrae. The 
penumbral distributions follow radial directions coming out from 
the umbra. Moat flows are organized radially as the penumbra. 
The coordinates are expressed in arc sec. A peculiar region is 
white squared. 

from top to bottom around coordinates (20,28) to (28,11) 
respectively (see Figure 2 in De Pontieu et al. (2007)). 
The penumbral filaments extend from the umbra to the 
right and seem to bend, being forced to follow the di- 
rection of the neutral line as they approach it. The 
sheared configuration of this penumbra is arranged so 
that the penumbral filaments end up along directions 
parallel to the sunspot border. Moat flows are then 
not found beyond this sort of penumbral configuration as 
mentioned in § 3.1. A similar case was also found by Var- 
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Fig. 6. — Sunspot S6: Map of the horizontal velocities inside the 
moat with de-projcctcd magnitudes > 0.3 kms^ 1 . This sunspot 
exhibits penumbrae on the right hand and upper sides. The penum- 
bra on the right is mostly tangential to the sunspot border. In the 
upper penumbra we identified the moat flow as the prolongation 
of the pcnumbral filaments as expected. The coordinates are ex- 
pressed in arc sec. A peculiar region is white squared. 
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Fig. 7. — Sunspot S7: Map of the horizontal velocities surround- 
ing the sunspot with de-projected magnitudes > 0.3 kms -1 . The 
sunspot includes a well developed penumbra but also an empty 
region with no penumbra. Looking at the large flows, they com- 
pletely dissapear in the emptied areas with no prenumbra. The 
white contours outline the moat regions. The coordinates are ex- 
pressed in arc sec. 

gas Dommguez et al. (2007) in a complex (^-configuration 
active region. For that active region the authors found a 
strong sheared neutral line crossing a penumbral border 
where a moat flow was expected to follow the penum- 
bral filaments direction but actually not detected. More 
observations of complex active regions with neutral lines 
present in the vicinity of penumbrae are needed to firmly 
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Fig. 8. — From top to bottom: Close-ups for peculiar regions in 
sunspots SI, S4, S6 and S5 respectively. The upper three regions 
show a penumbra extending tangential to the sunspot border. The 
lower pannel shows a radial penumbra. The coordinates are ex- 
pressed in arc sec. 
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establish the relation between the absence of moat flows 
and magnetic neutral lines. 

3.3. Statistics of velocity fields in moats vs. quiet 
granulation 

In this section the statistics of time-averaged horizontal 
velocity fields in moats is performed. For the sake of com- 
parison, the velocity statistics in quiet granulation areas 
has also been computed. To that aim, boxes (~ 9 x 9 
arcsec 2 ) in regions of less-magnetized (quiet) granulation 
and far from the moat flows have been manually selected 
in every FOV. Table 3 summarizes the statistical proper- 
ties of the de-projected magnitudes of horizontal veloc- 
ity fields within both moat masks and selected granular 
boxes. Here we again remark that because of the possi- 
ble underestimation of the displacements by using LCT, 
we intend to perform not an absolute but a comparative 
statistical study between moat regions, where we find ra- 
dially organized outflows, and quiet granulation regions 
far from the sunspots. 

For each sunspot of our sample three velocity magni- 
tude maps result from averaging over different time pe- 
riods, namely: 5 min, 10 min and the whole duration 
of the corresponding time series . Table 3 shows that 
in most cases the mean velocity magnitude in moats is 
greater than in granulation. In a few cases and only for 
averages in short time periods (5, 10 min), we find similar 
values in both or even slightly lower velocities in moats. 
However, when averaging over a long time period, the 
difference of mean velocities in moats and granulation is 
in all cases positive and more conspicuous than for short 
time period averages (see column 10 in the Table). A 
similar behavior is obtained for the rms parameter (see 
column 9 in the Table). The maximum velocity values 
are also systematically larger in moats than in granu- 
lation. The described statistical behaviour is expected 
since for short time averaging periods, of the order of 
the granulation lifetime (5-16 min; see Hirzberger et al. 
(1999) and references there in), the proper motions of 
the granulation structures compete in magnitude with 
the velocity of large-scale flows. However, averaging ve- 
locity components over long time periods results in lo- 
cal velocity cancellations in short-lived structures while 
steady motions at large spatial scales prevail. So moats 
can be considered as a long-term and large-scale phe- 
nomenon where the mean velocity exceeds that of quiet 
granulation by about 30%. 

The significant fluctuations in A (rms) and A (mean) 
(columns 9 and 10 in Table 3) for the velocity fields av- 
eraged over the whole time series, deserve a comment and 
possibly future improved measurements. They could be 
ascribed to morphological differences in the various cases 
considered or they could be related to the particular evo- 
lutionary state or magnetic field strength in the various 
sunspots of our sample. Furthermore, the accuracy in 
the measurement of 9 and <f>' in equations (1) and (2) is 
crucial for the determination of v. Thus, two sources of 

7 Actually, when averaging over 5 min (or 10 min) we obtain 
more than one map since we then compute maps every 5 min (or 
10 min) up to complete the total duration of the series (i.e. when 
having 40 min of duration for a time series, we would finally obtain 
8 maps of 5 min or 4 maps of 10 min). This enables us to improve 
the statistics of the velocity magnitudes and reduce the noise in 
our calculations. 



inaccuracy are: a) the assumption of constant heliocen- 
tric angle 9 all along the sunspot area (this has important 
impact in regions far from the solar disc centre); and b) 
the determination in our FOV of the direction pointing 
to the solar disc centre which defines the X-axis of the 
Observing Coordinate System, i.e. the angular origin for 
</>'• 

Figure 10 shows the histograms of de-projected veloc- 
ity magnitudes for averages over the whole time series for 
sunspots SI to S7. Thick lines correspond to velocities 
inside the moats and the thin ones to velocities in quiet 
granulation boxes. In all erases we find the same gen- 
eral trend in both histograms. The thick histograms are 
globally shifted toward the right with respect to the thin 
ones. The left wings of the histograms for granulation lie 
systematically above those of the histograms for moats, 
and for SI, S2, S3, S4 and S6, both histograms nearly 
coincide at the most left end. This means that within 
the moats standard velocity values for quiet granulation 
are still present but with a lower weight. Such is the case 
for the fragments of exploding granules that additionally 
are swept by large-scale moat flows. At some point from 
left to right, both histograms cross each other so that 
the right wing of the moat histogram surpasses the cor- 
responding quiet granulation wing and extends to larger 
velocity values (> 0.6 kms -1 ). This confirms the pre- 
dominance of large velocities in moats. The intersection 
point of both histograms corresponds to ~ 0.3 kms -1 in 
almost all cases (a bit higher in SI, S2 and S4). Interest- 
ingly, this value coincides with the threshold empirically 
selected a priori ( § 3.1) to clearly distinguish the fron- 
tiers of the moats, which supports the goodness of the 
limit- value (0.3 kms -1 ) chosen to detect moats in our 
data set. 

4. DISCUSSION 

Seven time series of sunspots were restored from instru- 
mental and atmospheric aberrations using MFBD and 
MOMFBD techniques. The high quality of the images 
allowed us to study proper motions of granules outside 
the sunspots and measure their time-averaged velocities. 

We have extended the study by Vargas Dominguez et 
al. (2007) to a larger sample of active regions and sys- 
tematically confirmed their findings: a) moat flows are 
oriented following the direction of the penumbral fila- 
ments; b) in granulation regions found adjacent to an 
irregular penumbral side parallel to the penumbral fila- 
ments moats are absent, or in other words, moats do not 
develop in the direction transverse to the penumbral fila- 
ments. Note that if the moat flows were originated by the 
blockage of the heat flux from below by the penumbra, 
one would expect moat flows directed along, but also, 
transverse to the direction of the penumbral filaments; 
c) umbral core sides with no penumbra do not display 
moat flows. Moreover we include in our sample a case in 
which a neutral line extends along a penumbral border 
where we would expect a moat flow continuation. For 
this sunspot we do not find any moat flow following the 
direction of the penumbral filaments after crossing the 
penumbral border where we see a change in magnetic 
polarity. The same result is found by Vargas Dominguez 
et al. (2007) in a penumbral portion of a complex active 
region crossed by a strong sheared neutral line. 

All these results indicate a likely connection between 
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TABLE 3 

Statistics of de-projected horizontal velocity magnitudes [m s -1 ] 



Sunspot 


Duration 




Moat 




Granulation 


A (rms) 1 


A(mean) 1 




(min) 


max 


rms 


mean 


max 


rms 


mean 


70 


0/ 
70 




5 


1048 


174 


333 


947 


157 


315 


10.8 


5.7 


C!l 
ol 


10 


947 


162 


315 


842 


157 


314 


3.2 


0.3 




53 


639 


137 


301 


596 


123 


250 


11.4 


20.4 




5 


1187 


183 


329 


898 


162 


307 


13.0 


7.1 




10 


1278 


182 


327 


862 


155 


300 


17.4 


9.0 




A *J 

4/ 


921 


150 


303 


678 


115 


249 


30.4 


21.7 




5 


1036 


146 


290 


870 


147 


282 


-0.7 


2.8 


CI"} 


10 


868 


146 


281 


671 


137 


289 


6.6 


-2.8 




A *J 

4/ 


689 


126 


265 


451 


101 


220 


24.8 


20.5 




5 


1060 


198 


411 


863 


164 


362 


20.7 


13.5 


o^t 


10 


954 


180 


379 


/ OU 




^AR 
o^o 


16.9 


8.9 




45 


834 


179 


365 


655 


118 


262 


51.7 


39.3 




5 


999 


178 


355 


968 


172 


307 


3.5 


15.6 


S5 


10 


972 


187 


372 


1017 


159 


291 


17.6 


27.8 




40 


802 


141 


355 


545 


107 


220 


31.8 


61.4 




5 


699 


129 


257 


683 


130 


278 


-0.8 


-7.6 


S6 


10 


673 


131 


248 


723 


136 


277 


-3.7 


-10.5 




40 


622 


125 


243 


526 


104 


210 


20.2 


15.7 




5 


872 


158 


300 


759 


142 


226 


11.3 


32.7 


S7 


10 


785 


159 


320 


828 


149 


227 


6.7 


41.0 




40 


692 


131 


285 


593 


121 


196 


8.3 


45.4 



tA(rms) and A(mean) stand for increments of the rms and mean 
values in moats with respect to quiet granulation. 

the moat flows and flows aligned with pcnumbral fila- 
ments. In a recent work, Cabrera Solana et al. (2006) 
suggest that the Evershed clouds inside penumbrae prop- 
agate to the surrounding moat and then become MMFs 
after crossing the sunspot border. The MMFs displace- 
ments trace very well-defined paths that can actually be 
clearly seen when averaging magnetograms in time (Sainz 
Dalda & Martinez Pillet 2005). Some of these MMFs are 
seen to start inside the penumbra. 

In this paper we also complement the study by Var- 
gas Dommguez et al. (2007) with a statistical analysis 
describing the differences between velocity fields in moat 
flows and in less-magnetized solar granulation nearby. 
In contrast to the granulation, moat flows are well- 
organized, steady and large-scale motions. For averages 
made over more than 40 minutes, the mean velocity in 
moats (0.3 kms -1 ) exceed that of quiet granulation (~ 
0.23 kms -1 ) by about 30%, although we obtain a con- 
siderable dispersion in the results. Also the rms of the 
velocity magnitude is greater in moats by a similar per- 
centage. 

The histograms of velocity magnitudes in the moats 
are broader than those in granulation. The histograms of 
granulation show conspicuous maxima, most of them in 
a range from 0.2 to 0.3 kms -1 , whereas the histograms 
of moats present a flatter top. Systematically, the low 
velocity wing of the granulation histogram lies above that 
corresponding to the moat. At some point about 0.3- 
0.4 kms -1 , both histograms cross each other and the 
right wing of the moat histogram extends beyond that of 
the granulation to larger velocity values (> 0.6 kms -1 ). 

We have studied one case of a sunspot penumbra dis- 
playing a neutral line all along a sector of the penumbra 
(Figure 9). This neutral line is detected at the penumbral 
boundary but also emphasized by a large opposite polar- 
ity concentration nearby. Interestingly, this penumbral 
sector shows no moat flow. We interpret this evidence 



as an indication of the Evershed flow being forced to go 
into deeper subphotospheric layers at a faster pace than 
what is normally though to occur in penumbral regions 
not associated to neutral lines (Westendorp Plaza et al. 
1997). A similar case, but in a sheared neutral line of a 
5-spot, was found by Vargas Dommguez et al. (2007). 

Though there is increasing evidence linking the moat 
flows and the Evershed flow along the penumbral fila- 
ments, the debate regarding the existence of a moat flow 
around umbral cores and individual pores is still undergo- 
ing. In a recent work, Deng et al. (2007) , found that the 
dividing line between radial inward and outward proper 
motions in the inner and outer penumbra, respectively, 
survived the decay phase, suggesting that the moat flow 
is still detectable after the penumbra disappeared. 

Previous works (Sobotka et al. 1999; Roudicr, Bonet 
& Sobotka 2002), have measured horizontal proper mo- 
tions in and around pores and have observed a ringlikc 
structure of positive divergence ("rosettas") around the 
pores, which is related to a continuous activity of explod- 
ing granules. Roudier, Bonet & Sobotka (2002) identi- 
fied a very clear inflow around pores which corresponds to 
the penetration of small granules and granular fragments 
from the photosphere into the pores, pushed by granu- 
lar motions originated in the divergence centres around 
them. They conclude that the motions at the periphery 
of the pore are substantially and continuously influenced 
by the external plasma flows deposited by the exploding 
granules. We interpret the dividing line between radial 
inward and outward motions, found by Deng et al. (2007) 
outside the residual pore, as corresponding to the centres 
of divergence of the exploding granules around the pore. 
The outward motions these authors described, which are 
not in the immediate surroundings of the pore but sep- 
arated by the annular inward motion would then corre- 
spond to the flows coming out from the regular mesh of 
divergence centers around the pore. 
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Fig. 9. — Sunspot S3: Upper panel: Map of the horizontal veloci- 
ties with deprojected magnitudes > 0.3 kms" 1 . The sunspot in the 
figure is surrounded by a well-developed penumbra. Large outflows 
are not found in part of the right-hand side of the spot. When look- 
ing at the magnctogram for the same field-of-view (lower panel) we 
found a neutral line which is suspiciously acting as a blocking agent 
for the flow coming out from the spot. The black box is a common 
region for both figures. The coordinates are expressed in arc sec. 
The black bar at (0,0) in the flow map represents 1.5 kms - 1 for 
the projected velocities. 

The important questions related to the flows inside and 
outside sunspots are yet being studied. In the present 
work we contribute with a new sample of sunspots ob- 
served between 2003 and 2006 in six different observing 
campaigns. The data sets have been selected on basis 
of the seeing quality (sharpness, homogeneity, duration) 
and on the availability of suitable targets: the presence of 
spots with some form of irregular penumbra. In all of our 
samples we follow the evolution of the sunspots for more 
than 40 min. Although this only represents a snapshot in 
the evolution of the sunspots through all their emerging 



and decaying processes, our sample includes sunspots in 
different evolutionary stages and penumbral configura- 
tions. 

New facilities such as the recently launched HINODE 
satellite can provide long time series of active regions 
with a constant image quality and enough spatial resolu- 
tion to provide firm confirmation of the evidences found 
in this paper. The addition of simultaneous Dopplcr 
and magnetogram data to the continuum intensity (or 
G-band) data sets will enhance our understanding of the 
link between the Evershed and the moat flows. Need- 
less to say, the study of those stages were penumbrac 
are just formed or destroyed becomes of extreme impor- 
tance to validate our findings. Similarly, the results that 
can be expected in the coming years from local helioseis- 
mology, describing the flow patterns in the deeper layers 
near sunspots, will prove crucial for the establishment of 
a clear link between these two well known flow patterns 
that has so far been not appreciated. 
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the Spanish Ministerio de Education y Ciencia through 
project ESP2003-07735-C04 and financial support by the 
European Commission through the SOLAIRE Network 
(MTRN-CT-2006-035484) are gratefuly acknowledged. 
This research was supported through grants 146467/420 
and 159137/V30 of the Research Council of Norway. 




Fig. 10. — Histograms of the velocity magnitudes in moats (thick line) and quiet granulation (thin line), averaging over the whole time 
series (more than 40 minutes in every case) for sunspots SI to S7. The vertical and horizontal axis represent the percentages and the 
velocity magnitudes [kms" 1 ] respectively. 
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